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Abstract 

We present a new perturbation theory for quantum mechanical energy eigenstates when the 
potential equals the sum of two localized, but not necessarily weak potentials V\{r) and V^(r), 
with the distance L between the respective centers of the two taken to be quite large. It is assumed 
that complete eigenfunctions of the local Hamiltonians (i.e., in the presence of V\{r) or V2{r) only) 
are available as inputs to our perturbation theory. If the two local Hamiltonians have degenerate 
bound-state energy levels, a systematic extension of the molecular orbital theory (or the tight- 
binding approximation) follows from our formalism. Our approach can be viewed as a systematic 
adaptation of the multiple scattering theory to the problem of bound states. 
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I. INTRODUCTION 



In one-particle quantum mechanics, consider the energy eigenvalue problem 

£|V> = £|V), (i) 

when the Hamiltonian of the system has two separate potential contributions, viz., 

H = ^-p 2 +Vx + %. (2) 
2m 

(For simplicity, we will present our discussion within a one-dimensional context). Then one 
may also consider the related eigenvalue problems 

£ihM = e|Vi>, (3) 
H 2 \ip 2 ) = u\ip 2 ), 



where 



H 1 = 7 ^-p 2 + V 1 , H 2 = ^-p 2 + V 2 . (4) 
2m 2m 



Generally speaking, there will be no simple connection between the eigenvalue problem ([!]) 
and those in (|3]) (other than an inequality-type relation for the ground state energy). But, if 
V± and V 2 correspond to some localized, but not necessarily weak, potentials with the centers 
at x = and x = L, respectively and the separation distance L is relatively large, one might 
hope that the solutions to the eigenvalue problems in ([|) (involving 'local Hamiltonians' 
Hi and H 2 ) be useful for generating good approximate solutions to the initial problem ([[]). 
Indeed, this view forms the basis of the so-called molecular orbital theory or the tight-binding 
approximation [1,2], in which one diagonalizes the full Hamiltonian H within the truncated 
vector space given by a linear combination of atomic orbitals (consisting of a few low-lying 
eigenstates of the local Hamiltonians). In the context of Born-Oppenheimer approximation 
where the local potential centers are not really fixed, this kind of energy eigenvalue problem 
is of particular importance since it can account for an effective binding force between the 
potential-producing objects. 

The tight-binding method or its variants will be useful when given local Hamiltonians 
allow some deeply-bound orbitals which are separated from other local eigenstates by rel- 
atively large energy gap. By its very nature, however, a reliable theoretical error estimate 
for the scheme (especially when the parameters in the given problem are not quite in the 
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limiting range for the method) is difficult to make. Also, if one of the local potentials, say, 
V 2 , happens to be strictly repulsive (and so no atomic orbital associated with H 2 ) while 
Hi allows some bound states, this method is unable to give any useful information on the 
effect of the potential V 2 on the low-lying eigenstates of the full Hamiltonian (0). There is 
a related question within the usual tight-binding approximation, that is, on the role of the 
continuum states in the scheme. Clearly, it is desirable to have a systematic approximation 
scheme which goes beyond the simplest tight-binding approach. [Recently, Barton et al.[3] 
discussed the effect of a distant impenetrable wall on quantum mechanical energy levels; 
but their approach is tuned to the change of the boundary condition, and therefore does not 
apply to more generic case involving two well-separated potentials.] 

In this paper we develop a new stationary perturbation theory which can be used to 
study the eigenvalue problem with a two-centered Hamiltonian. (For a Hamiltonian with 
more than two centers a simple extension of our method should be useful.) While there ex- 
ists a systematic theory dealing with scattering by a multi-centered potential (see Ref. [4] for 
instance), we are not aware of such development which can be used to study the correspond- 
ing bound-state problem in a well-controlled manner. In our approach to the eigenvalue 
problem ([[]), it will be assumed that the eigenvalue problems with the local Hamiltonians 
can be solved explicitly, and so we have at our disposal a complete orthonormal set {\n)} 
based on eigenstates of Hi and another complete orthonormal set {\n)} based on eigenstates 
of H 2 . [The knowledge of the Green's operators associated with the local Hamiltonians may 
be assumed instead.] We wish to exploit this over-complete set of basis, which include con- 
tinuum states, in constructing the bound states of the total Hamiltonian H. The result is 
a perturbation series in which the expansion parameter is a quantity approaching zero as 
the separation between the local potentials becomes large. [In fact, for strongly localized 
local potentials, we have an expansion parameter of order e~ aL (a: constant)]. It can be 
viewed as an expansion in the wave-function stretching factor, that comes with for every V 2 
(Vi) acting on a specific bound state of Hi (H 2 ). This small factor is a direct measure on 
how much influence one local potential feels from the bound states associated with the other 
local potential. We also remark that the general philosophy of our formalism is similar to 
that of the multiple scattering theory but the very nature of the bound-state eigenvalue 
problems necessitates somewhat different developments. 

It should be noted that the standard time-independent perturbation theory is generally 
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unreliable for our problem. To see that, it suffices to consider the simple situation where 
one has the bound state energy levels of Hi influenced by a strictly positive, well-localized 
potential V 2 (x) at a large distance L (from the center of the potential Vi(x)). The strength 
of V 2 may not be small, however. If this case can be studied by the usual perturbation 
theory (that is, by treating V 2 as a perturbation to the unperturbed Hamiltonian Hi), the 
state \k) satisfying Hi\k) = Ek\k), i.e., that with the unperturbed energy = at (< 0), 
would acquire the first- and second-order energy shifts 

E W = (k\V 2 \k), (5a) 
E {2) = (k\V2\n)(n\V 2 \k) ^ 

assuming for simplicity no degeneracy for the unperturbed states. According to (|5a|), E^ 
would be of order e _E "v/ 2m ' £fc ' L . This is nothing but the product of two wave- function 
stretching factors, as appropriate to the matrix element of V 2 in a specific bound state of 
Hi. The fact is that, according to £ (2) would also be e>(e _ V^"^ L ) due to the 

continuum contribution in the intermediate-state sum. This implies that, depending on the 
strength of V 2 , the second order shift might be as big as the first-order shift. In an 
analogous manner, it is not difficult to see that the contributions from the continuum states 
make the r-th order shift E^ assume the same order of magnitude as E^\ Hence this is 
not a valid expansion, and we have to devise a more elaborate scheme to solve our problem. 

This paper is organized as follows. In Sec. 2, we will concentrate on setting up a reliable 
perturbation theory with two spatially well-separated potentials in the nondegenerate case. 
In this discussion we will suppose (mainly to have mathematics under control) that the 
potentials Vi and V 2 are sufficiently well localized; but, we expect that most of our formulas, 
with suitable adjustments if necessary, remain useful even if these potentials are localized 
only by some (not too small) powers in the distance from the respective potential centers. 
Our method is exhibited explicitly for local Hamiltonians involving 5— function potentials. 
Section 3 is devoted to the extension of this method to the case where the local Hamiltonians 
Hi and H 2 have (almost-)degenerate energy levels. Here one sees explicitly that, for a reliable 
perturbation series, a separate treatment in the subspace of degenerate local bound states 
becomes necessary. The resulting theory is a generalization of the molecular orbital theory 
that allows one to systematically study higher order corrections, and as such it should have 
some practical value as well. Section 4 contains concluding remarks. In the Appendix we 
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present our argument behind the order estimates for various contributions appearing in our 
perturbation theory (together with some analysis for the example problem). 

II. NONDEGENERATE PERTURBATION THEORY 

Our goal is to obtain approximate eigenstates of the Hamiltonian H = j^p 2 +Vi+V2, when 
complete solutions to the eigenvalue problems with the local Hamiltonians H±(= + Vi) 
and H 2 (= ^mP 2 + ^2) are known. In the Hilbert space V of the system we have with us two 
complete orthonormal sets — the set {\n)} based on (discrete and continuous) eigenstates of 
H\ and the set {\n}} based on eigenstates of H 2 . Let \k) be a given specific nondegenerate 
bound state of Hi, with eigenvalue e k . We further assume in this section that no eigenstate 
of H 2 has the eigenvalue equal or very close to e k . Then, if the distance L between the centers 
of two local potentials Vx(x) and V 2 (x) is large enough, we expect that the full Hamiltonian 
H admit an energy eigenstate | </>/-) which should coincide with \k) in the limit L — > 00 (i.e., 
as V2(x) is sent away to the very remote). This should be the case irrespectively of the 
relative magnitude of the two local potentials. Thus, for large L, we may write the solution 
to the eigenvalue equation 

H\<p k ) = E k \<p k ) , {H = H x + %) (6) 

as 

E k = e k + 5E k , \<f> k ) = \k) + \S<f> k ). (7) 
Here the small corrections 5E k , \5<p k ) should satisfy the equation 

(s k - H + 5E k )\5<p k ) = (V 2 - 5E k )\k), (8) 

which is still exact. 

Let us now study the implication of in detail. First of all, as in the ordinary stationary 
perturbation theory, (@) does not determine \84> k ) uniquely [5]: if \5<pk) is a solution of (f|), so 
is \54>k)' — i^{|^0fc) — P\k)} for arbitrary constant /3. As a result, {k\S(f) k ) may be chosen 
as one wishes and the particularly convenient, at least in the ordinary perturbation theory, 
is the choice 

(k\5(f> k ) = 0, (9) 
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i.e., define \6<j>k) in the subspace the orthogonal complement of \k) in V. With this choice 
and multiplying both sides of (^|) by (k\ on the left, one obtains 

5E k = (k\V 2 \k) + (k\V 2 \5<p k ) (= (k\V 2 \cj> k )) . (10) 

At the same time, one may replace (|8|) by 

Q k {ek-H + 5E k )\5<f> k )=Q k V 2 \k), (11) 

where Q k = 1 — Note that, without the knowledge on \5(f> k ), the formula (|10D is not 

informative by itself. To have determined, one might write (as in the conventional 

perturbation theory) \8<p k ) = J2n{^k) \n)(n\S<pk) and determine (n\S<p k ) with the help of the 
equations resulting from multiplying (|TT|) by (n\ on the left. But, as was explained in the 
introduction, this usual procedure does not lead to a useful perturbation series. (See also 
discussions further below.) 

At this point, recall that, in association with the second local Hamiltonian H 2 , we have 
another complete set {\n)} where H 2 \n) = Un\n). We shall utilize them with (ffl) in a suitable 
manner. (See (|16j) below.) Here it is convenient to recast (||) and ( |TT| ) as 

6\S<f> k ) = V 2 \k)-SE k \k), (12) 

QkO\5<f> k ) = QkV^k), (13) 

introducing the operator 

6 = e k -H + 5E k = e k -H 2 -(V 1 -5E k ). (14) 

But we are not going to use the condition (Q) — it is not convenient for our development. 
[Note that (11) holds good without assuming this condition]. On \5<f) k ) we only demand 
that it should be small, i.e., suppressed by at least one wave-function stretching factor 
(accompanying, say, a term like V 2 \k)). If we multiply (|T2"D by (k\ on the left without 
imposing (|9|), we obtain 

Ap (k\V 2 \k) + (k\V 2 \6<p k ) 

From this formula, we may conclude that 5E k contains at least two wave-function stretching 
factors. Now note that, when G 2 denotes the Green's operator associated with the second 
Hamiltonian H 2 

= (16) 



the operator O satisfies the relation 

1 = OG 2 + (V a - 6E k )G 2 , ■ (17) 

Hence the right hand side of (|T3| ) may be written as Q k {OG 2 + {V\ — 8E k )G 2 }V 2 \k) , and 
then, by rearranging, we obtain 

Q k O {\5<j) k ) - G 2 V 2 \k)} = Q k V x G 2 V 2 \k) - 8E k Q k G 2 V 2 \k). (18) 

This equation is the crucial one for our perturbation scheme. 

We wish to solve flT8|) order by order, with the order in our case determined by the number 
of the wave-function stretching factors involved. In its left hand side we have the operator O 
acting on a vector yet to be found, \w) = \5(p k ) —G 2 V 2 \k). Here it is important to note that, 
if |w) does not contain a component proportional to \k), |w) and 0\w) would be of the same 
order due to the assumed nondegenerate nature of H. As for the component proportional 
to \k) from |w), on the other hand, the situation is not the same: if O acts on that piece, 
the resulting vector will have the order increased by at least one wave-function stretching 
factor. This follows from 

6\k) = {e k -H 1 + 8E k -V 2 )\k) = (8E k -V 2 )\k). (19) 

Based on this observation, the following conclusion should be immediate: in (0), the vector 
\w) (= \5(j) k ) — G 2 V 2 \k)) appearing in its left hand side is necessarily of the same order 
as the expressions in its right hand side, under the proviso that this restriction on the 
order does not apply to the term proportional to \k). We here make another important 
observation: the expressions we have in the right hand side of (|TH|) are in fact of higher 
order than that of G 2 V 2 \k). For its justification, see the Appendix. Hence, setting \w) ~ 0, 
i.e., \5(f> k ) — G 2 V 2 \k) solves ( |18|) to the leading order. We may thus write 

\6<t> k ) m =G 2 V 2 \k) (20) 

and, using this with (0), the following formula for the energy shift results: 

5E ( k 1] = (k\V 2 \k) + (k\V 2 G 2 V 2 \k) . (21) 



[Note that (k\5(j) k ) in the denominator of (15) can be ignored in the leading-order approxi- 
mation.! 



With the expression ( |i~6"D for G 2 inserted into (p0|) and fl2"l|) , our formulas giving the 
leading-order correction to the energy eigenstate may be recast as the ones involving the 
sum over the basis set 

W» = EI«P- (22) 
5 E M = {k \v 2 \ k) + E . (23) 

Note that, in our procedure, no explicit condition (like that in (|9])) has been used to dispense 
with the ambiguity concerning the |fc)-component of \5(f) k ). Instead, we have decided to 
choose the simplest available expression for \S(f> k ), as suggested by the order- by-order analysis 
of the relevant equation for our perturbative development. In view of (^), one may well 
say that our choice in fact corresponds to 

(k\5(j) k ) = (k\G 2 V 2 \k) + {higher order). (24) 

The energy eigenstate we obtain is not properly normalized in general. 

Suppose one attacked the above problem with the help of the conventional perturbation 
theory, regarding V 2 as a would-be perturbation. Then, instead of ([18]), one would work 
with the equation 

Q k 6{\5<p k ) - G[V 2 \k)) = QkVzG'&lk) - 6E k Q k G[V 2 \k) . (25) 

where G[ = J2 n (^k) 1S the Green's operator associated with the local Hamiltonian 

H\ (but defined in the orthogonal complement V k ). To obtain (^), one may utilize the 
equation (|35"D given below with (|13|). The lowest-order approximation in the conventional 
perturbation theory is tantamount to identifying \5<f) k ) with the term G^V^k). But, in our 
case, this is not a good approximation (unless the strength of V 2 itself is very weak), since 
the first term in the right hand side of ( Po]) can generate a comparable contribution. (See the 
related discussion in the introduction). Note that we had a different situation with flT8| ) - 
the expressions in its right hand side were of higher order (i.e., involved more wave-function 
stretching factors)! 

To be convinced of the validity of our leading-order approximations in (^) and (p3|), let 
us consider a simple example consisting of a pair of 5-function potentials, i.e., 

V x {x) = - 7 i5(x) , V 2 (x) = - l2 5(x - L) (26) 
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with 71 > 72 > 0. Then we know that each local Hamiltonian admits one bound state. If 
|1) (|2)) denotes the bound state of H\ = ^p 2 + V\ (of H 2 = -^p 2 + V2), the corresponding 
(normalized) wave-function and energy eigenvalue read 



M*) = M 1) = e-^W , (with e, = (27) 



6,(«) S (x|2> = ^ e-^*-« , (with u 2 = (28) 

respectively. Now, if the distance between the two local potentials, L, is large (and the 
value of U2 differs from that of E\ significantly), we expect that the full Hamiltonian H = 
2^j3 2 + V1 + V2 allow two bound states \<f>i) and \<f>2), which are approximately equal to |1) 
and 1 2), respectively. For this example one can of course find the exact bound state energies 
by a direct analysis of the corresponding Schrodinger equation. Explicitly, for the state \4>i), 

2 

its energy E = —^r is determined by the equation 

rf - (71 + 72)?7 + 7i72(l - e- 2 ^ r >) = 0, (29) 
and therefore, for large L, one has 



mjf f 272 _ 2 nm 



The above result can also be obtained by using our formula (P3"|). For such check, we 
need a complete basis {\n)} consisting of the energy eigenstates of H 2 - in position space, 
the desired complete set contains, aside from the bound state £0(^)5 two distinct classes of 
continuum states (both corresponding to energy u q = -o£~) 



/ l-fTh 1 r 1 rl -1/^72x1 

{x\q v ') = —= cos|g|x - L\ + tan (— T )\, 
Y vr ~ qh 

(x\q^) = —= sin[g(x — L)]. 

'71 



(31) 



Then, by straightforward calculations using these eigenfunctions, we find 



<1|V 2 |2)<2|VJ|1) _ 4(5; 



,2' 



si-u 2 (2i)2_i v 2^ 2 



m71 * e" 2 ^, (32b) 



..... r, (i|^|g (1) )(g (1) |v r 2|i) + (i|v r 2|g (2) )(g (2) |v r 2 |i) 2(g) (m 7 2 \ 

(m) l dq — q =-(^ri{-w) e R - (32c) 
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By summing these contributions, we thus obtain the result 8E^' = — (-TrS-le - ~& 

J a i 71—72 v 2h 1 



which is in agreement with (|30|). One may also calculate the first order eigenfunction cor- 



rection with the help of our formula (|22|). After some straightforward calculations, we then 
find 

m7l f -SXLM 72 -^L -^Lb-ill ..,->, 

J e 1 H e s 2 ^ e R^ 1 >. (33) 

71-72 I 



This is the correct result, for the exact eigenfunction in the limit e ~rF ^ (with 77 
= 7i + 0(e~ 2 ~R 2_i ) from ( p9|) ) can be approximated by 



^™7i | e -^N + _^_ e -^V^I- £ l| • (34) 

Expressions for higher order terms of our perturbation theory can be found also. Here we 
shall concentrate on identifying the second order terms, since even higher order terms can 
be found by a rather obvious extension of this procedure. For the purpose, we had better 
rewrite the contributions in the right hand side of flT5p appropriately. As regards the first 
term, we may here use (instead of (|T7|)) the identity 



Q k = OG[ + (V 2 - 8E k )G\ , f G[ = £ -MM. ] (35) 

to have it rewritten as 

Q k ViG 2 V 2 \k) = QkOG'^GzV^k) + Q k (V 2 - 8E k )G' x ViG 2 V 2 \k) . (36) 
On the other hand, with the second term, the relation (|i~7|) can be used to write it as 

- 8E k Q k G 2 V 2 \k) = -8E k Q k 6{G 2 fV 2 \k) - 8E h Q h (V-y - 5E k )(G 2 ) 2 V 2 \k). (37) 



Using these forms with (|T8| ) leads, after some rearrangements, to the following equation: 



Q k O {\5<p k ) - \6</> k )W - G[ViG 2 V 2 \k) + 5E k {G 2 fV 2 \k)} = 

Q k (V 2 - 8E k )G[ViG 2 V 2 \k) - SEkQ^ - 5E k )(G 2 ) 2 V 2 \k). 

In the Appendix the expressions in the right hand side of (|38|) will be shown to be of 
higher order than the terms appearing inside the curly bracket in the left hand side of the 
same equation. Note that, for this behavior, it is crucial to have Q k ViG 2 V 2 \k) rewritten as 
in (|36| ) and not by the form 

QkViG 2 V 2 \k) = Q k OG 2 ViG 2 V 2 \k) + Q k {Vi - 8E h )GAG 2 V 2 \k) , (39) 
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as the use of ([[7]) would result in. [Here the term involving ViG 2 V\ is dangerous, when con- 
tinuum contributions are considered.] On the other hand, as for the term — SE^ QkG 2 V 2 \k), 
it is allowed to have (|37|) replaced by another relation obtained with the use of fl3~5[) — but, 
using (|37D (and hence the equation (j38|) ) leads to a simpler perturbation theory practically. 
Now, based on this order count for the terms appearing on both sides of ([38]), we are led 
to conclude that the expression inside the curly bracket may be set to zero in our present 
approximation. Note that this reasoning is entirely similar to what we used with (Jl8|). As 
a result, it is found that our second order approximation to \5(pk) can be identified with 

\Hk) {2) = G'&G^k) - 5El 1) G 2 G 2 V 2 \k). (40) 

By using this expression with (|H>D, one can obtain the corresponding formula for the second- 
order energy shift also. 

For the third or higher order approximation, one may repeat the above procedure. Clearly, 
the approximation at desired order follows immediately once one has the appropriate gener- 
alization of the equation like ([18]) or fl3~8|) . As we have explained above, such generalization 
can always be found by using the identities (|i~7D and fl35|) in a judicious way with the corre- 
sponding equation one order lower. For a useful guideline here, see the Appendix. 

III. (ALMOST-) DEGENERATE PERTURBATION THEORY 

Our perturbation theory in the previous section was developed under the no degeneracy 
assumption; that is, for a given unperturbed state \k) (an eigenstate of Hi with eigenvalue 
£fc), no other eigenstate of Hi or H 2 has the corresponding eigenvalue equal or very close 
to Efc. In this section we will dispense with this restrictive assumption. The perturbation 
theory to be developed below is applicable to the case when (almost-) degeneracy, within the 
spectrum of Hi or between the spectra of the two local Hamiltonians Hi and H 2 , is present. 
This consideration is especially relevant since many physically interesting problems, which 
were treated traditionally by the molecular orbital theory, do come with such (almost-) 
degeneracy due to symmetry or by other reasons. 

First, we focus on the case when there are two (almost-) degenerate states \k) and \k) 
with Ek ~ Uf.. That is, each local Hamiltonian has a bound state of almost identical energy. 
[ This happens especially if H 2 is related to Hi by a simple spatial translation, i.e., V 2 (x) = 
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Vi(x-L). ] In this case we expect that the exact eigenstate \4>k) of the full Hamiltonian (0) 
have large overlap with both \k) and \k), in accordance with the philosophy of tight-binding 
approximation. So we may set up our perturbation theory by writing 

\<f> k ) = \k)+b\k) + \S<f> k ) (41) 

where \8<f>k) is supposed to be small, but the constant b can be a priori 0(1). [ Of course, 
if Ek and Uk were not very close to each other, b would become much smaller than 1. J 
Inserting the form ([|l]) into the Schrodinger equation (|6]) then yields 

6\5<p k ) = (^-SE^ + b^ + u^-Ek-SE^lk), (42) 

where 8Ek = E k — 6k, and O is the operator introduced in (14). Here the unknowns are 



5E k , b and \S(pk), and ( 42|) contains all the conditions required of them. 

If we multiply both sides of ( (42|) by (k\ or (k\ on the left, we obtain two relations which 
can be used to determine 5E k and the constant b, given the knowledge on \5<f) k ). Explicitly, 
we may write them as two different expressions for SE k , i.e., 

Ap a + bT + (k\V 2 \5(l> k ) 

5Ek = l + bA + (k\S^ k ) ' (43) 
b(u k -e k ) + T + b(3 + (u- k - s k )(k\5<p k ) + (fc|K|^ fc ) 

6Ek = bTA + WA) ' (44) 

where a, (3, T and A represent the matrix elements 

a=(k\V 2 \k), (5 = {k\V x \k), 

~ - - (45) 

T={k\V 2 \k), A=(k\k). 

[We have used the fact that, when \k) (\k)) is a nondegenerate eigenstate of Hi (H 2 ), it is 
possible to take a, (3, T and A to be real]. To fix the constant b (for given \S(j) k )), one can 
thus solve the quadratic equation obtained by equating the two expressions in the right hand 
sides of (|43|) and (|44f) . Then, how can one determine the eigenfunction correction \S(j) k )7 As 
in the nondegenerate case considered in Sec. 2, an appropriate perturbation theory for \5(pk) 
may be set up by considering the restriction imposed by (^2|) on its components belonging 
to the space orthogonal to \k) or \k). Again, in the corresponding development, we will not 
impose any specific condition on (k\5<p k ) (or, if one wishes, on (k\5<p k ) + b(k\5<f) k )); following 
the order-by-order analysis, it should suffice for us to choose \54>k) to be a simplest available 
expression that is consistent with the equation ([42]) . 
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When the separation distance between the local potentials is taken to be large, the leading 
approximation in our approach corresponds to the standard molecular orbital theory. This 
can be seen as follows. In ( f4~3"| ) and (|44|), b is order 1 while \8(p k ) is supposed to contain 
at least one wave-function stretching factor. Also, in the limit we are considering, all four 
matrix elements in fl4"5p should be quite small; a and f3 contain two wave-function stretching 
factors, and V and A one wave-function stretching factor each. In view of the potential 
V2 present in its definition, V may be estimated to be of order |u^|A. Moreover, from the 
assumed almost- degeneracy of the two states, it should be natural to assume that 

\u~ k -e k \ < \u- k \. (46) 

With T ~ |itjfc|A, this implies u k — e k <C T/A also. Then, in the leading approxiamtion, we 
may set \8(j) k )^ = (i.e., \(j>k) = \k) + b\k) to this order) and replace the right hand sides of 
([43|) and ( [44]) by bT and fc ~^ fc ' 1+r , respectively. From these, we conclude that 



6<°) = ± ^ (^p) 2 + 1 , 6EW = rfcW . (47) 

These are what one would expect with the original Hamiltonian replaced by the 2x2 matrix 
Hamiltonian (in the space spanned by two atomic orbitals \k) and \k)) 

( 



(48) 



E k + a r 
\ r u- k + (3 

and with a, (3 ignored because they contain two wave-function stretching factors while T has 
one. In particular, if u k —e k T, that is, if two energies are very close, the expressions in ( f4TD 
tend to the familiar values in the exactly degenerate case, b^ ~ ±1 (i.e., \<p k ) ~ \k) ± \k)) 
and 8E$ = ±r. It may also be of interest to look at the case u k — s k ^> T, that is, when 
the two energy values are not very close to each other (although they are almost degenerate 
in the sense of (f4"ED). Then, from the two values given for b^°\ only one of them — that 
with the behavior — > as T approaches zero — may be chosen since we are seeking for 
a solution that reduces to \k) in the absence of the potential Hence, with u k — e k ^> T, 
we find from (|47|) the values b^ = — - — and 8E1P = p2 , which are the results we can 
infer also on the basis of our formulas (|22D and (|23| ) (i.e., the lowest-order results in our 
nondegenerate formalism) . 



For higher order corrections, one should look for an iterative solution of (42), as we did 



the same with (|J) in the nondegenerate case. Here, for successive iteration, we will make 
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use of the relation (instead of (|I7D) 

Q k = OG' 2 + (Vi - 5E k )G' 2 , (49) 

where Q k = 1 — |fc)(A;|, and 

G' 2 = 4^4 = £ -MM- . (50) 

The Green's operator G[, satisfying (|35|) , will be useful as well. But we will here proceed 
somewhat differently from the nondegenerate case by not utilizing a suitably projected 



version of (42) in making iteration; for the present (almost-) degenerate case, manipulating 



directly with ([42]) is more convenient. Now note that, thanks to fl35|) and (§^), the terms 



V 2 \k) and Vi\k) in (42) can be rewritten as 

V 2 \k) = Q- k V 2 \k) + \k)(k\V 2 \k) 

= OG' 2 V 2 \k) + {V x - 5E k )G' 2 V 2 \k) +V\k), (51) 

Vi\k) = QkVi\k) + \k){k\V l \k) 

= dG^V^ + ^-SE^V^ + lT-iu-k-s^lk). (52) 



Using these in ( f42"D and then collecting all terms involving the operator O explicitly, we 
obtain the following equation: 



OllS^-G'Mlty-bG'Mk)) = V 1 G' 2 V 2 \k)-SE k G' 2 V 2 \k)+bV 2 G' 1 V 1 \k) 

-M^G^ilA;) - {5E k -br(l- ^)}\k) - {b5E k - T}\k). 

Based on (|3|), we will now show that the leading approximation for \S(f> k ) can be taken 



as 



= G' 2 V 2 \k) +bG' 1 V l \k). (54) 

First note that, as in the nondegenerate case, the first four terms in the right hand side of 
(|53|) can be shown to be of higher order than the expression G' 2 V 2 \k) + bG'^V-^k). On the 



other hand, the last two terms in the right hand side of (|53D are explicitly proportional to 
I A;) or \k); they are present because we are not working with a projected equation. Still, we 
observe that these terms are also smaller than the expression G' 2 V 2 \k) + bG[Vi\k), if the 
lowest order values for 5E k and b (in (47)) are used. [Here remember that u ^~^ < 1.] Then, 
based on these and our earlier observation as regards the effect of the operator O (in that 
case with (|i~8D), the identification (|5lD can be made. 
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The expression (54}) , with b replaced by b^°\ may in turn be used in (f43[) and ( |44"D to find 
the second order energy shift and the value . The results, to the appropriate order 

in the wave-function stretching factor (but without making an expansion with respect to 
another small factor ^7^ ), read 



"fc- £ fc 



h (i) f ^ (fc|ViGiVi|fc)-<fc|V 2 ^ 2 V 2 |fc) , . 

1 2r 2r M 



(55) 



= + a + (ife|y 2 Gf a y 2 |A;) - rA(M 0) ) 2 . (56) 
Especially, with = e fc , i.e., in exactly degenerate case, (|56"D reduces to 

^, = wM±|M^ + <5M±|MM)_ w ^ W y^), (57) 

as the definitions for a, /3, T and A in fl45|) are used. According to this formula, the second 
order energy shifts for the two split states become identical. An explicit check for the validity 
of (0) may be made for our 5-function example (see (p6|)) with 71 = 7 2 = 7. According 
to the direct calculation based on fl2§|), we have 5E ( k 2) = (2^ + l) e" 2 ^. We 

have verified that this very result is reproduced when various terms in Q57D are explicitly 
evaluated. Also, as in the nondegenerate further rearrangement of ( p3|) may be 

considered to obtain the expressions for the next order contributions. But, because of the 
complications involved and because their usefulness is rather limited, we will not consider 
such further higher order terms. 

It is possible to generalize the above discussion to the case when there are more than two 
degenerate states, that is, Aq eigenstates {(fc^) : \i = 1, ■ • • , N{\ of Hi with the given energy 
Ek and N 2 eigenstates {\kp) : p, = 1, • • • , N 2 } of H 2 , with the same energy uj, = Here 
we will concentrate on exactly degenerate case, not to make the problem too complicated. 
Now, for the exact eigenstates of the total Hamiltonian, we may write 

|0>=£ a A> + EM^> + i^>> (58) 

where a M and bp can be 0(1), but \5<p) is small. Inserting this form into the Schrodinger 
equation (f|), we obtain an equation similar to (42|), 



6\8<j>) = Y, a ^-5E k )\k^) + Y,h{V l -5E k )\k fi ). (59) 
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Then, from multiplying both sides of this equation by {k^\ and {kp\ from the left, we obtain 
the following conditions which may be used to determine 5E, a M and b^: 

5E(a^ + Y j A^b v + (k^\54>)) = E r > & p + E c W^ + (k„\V 2 \S((>), (60) 

V V V 

5E(b fi + E 4^ + (h\ s <t>)) = E r k^ + £ A*A + (h\Vi\64>), (61) 

V V V 

where we have defined 

<V = (kt,\V 2 \K), = (kplVtlh), ,„„. 

(62) 

T M p = (k^Vilkp) = (kp\V2\kp), A mP = (kftlkp), 

and f denotes the hermitian conjugate. [Note that T and A are Ni x N 2 matrices — not 
square matrices in general.] 

To determine the lowest order values 5E^\ affl and b^ , we note that (|60|) and floTf) , as 
only leading order terms are kept, imply the following equations: 

^ (1) < 0) = E r ^ 0) > (A* = l,"-,^i) 

SE^bf =ErU 0) - 0s = i,---,Ai). (63) 

V 

These can be regarded as a single eigenvector equation for an (Ni + N 2 )- vector (a^°\b^), 

^ (;»>)• < m > 

This is equivalent to the molecular orbital theory approximation [1,2] in which the full 
Hilbert space is truncated to the finite-dimensinal space spanned by Ni+N 2 atomic orbitals, 
i.e., {\k^)} and {|^)}. Assuming iVj > N 2 , the Ni + N 2 eigenvectors and corresponding 
eigenvalues may schematically be expressed by the forms 

(af^b^) = (u I „±v I ,), 5E? ) =±\ I , for 7 = 1,2,- --AT, 
(«j>S) = ( f/ ^ ) 0), 5^ = 0, for J = 2N 2 + 1,---N l +N 2 . 

The first 2A^ 2 eigenvectors are given by N 2 pairs of states, i.e., {E^mI^m) ^E^aI^m) > ^ = 
1, • • • , AT 2 }, with respective energy splits ±Aj. [Here, from studying (0), it can be shown 
that ^2 u* I(J ui^ = ^2 v^^vjfj, for each I — 1,2, ■ ■ ■ , N 2 , and so all uj, Vj may be taken to be 

unit vectors.] If the eigenvalue set {±A*} contains zero or the same value more than once, the 
degeneracy is not completely lifted and one may have to perform higher order analysis for 
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the effect (and associated true energy eigenvectors). There is no first-order energy shift for 
the remaining Ni—N 2 eigenstates, represented by Uj^k^) ; J = 2N 2 + 1, • • • , N\ + N 2 }. 
Hence, with N± > N 2 , there always remains some energy degeneracy which is not lifted by 
the lowest order consideration alone. 

To develop the corresponding higher-order perturbation theory, one should now take the 
expressions 

as relevant Green's functions and proceed in more or less the same manner as in our earlier 
consideration. Especially, with Ni = N 2 = N, we then find the results (as direct generaliza- 
tions of (0) and (|3)) 

|50 7 }« = $3«^2^2|M±E v ^i^il^>. ( 67 ) 
5E\ 2) = ^Y. u %{^ + (K\V2G' 2 V 2 \k u )}u Iv + ^ 

_ tEK a ^+ i '^%} 1 ( 6s ) 

where we have normalized Uj and v i to be unit vectors. 

With N\ > N 2 , we need to consider also the higher order terms to determine the above 
Ni—N 2 eigenvectors Ui^k^) ; J = 2N 2 + 1, • • • , ATi+A^} unambiguously and the possible 
energy splitting between them. To that end, one has to study the second order contributions 
from (pD|). Let us here assume for simplicity that the (iV\ — A^-dimensional space spanned 
by the states V] Ui^k^) represent the entire subspace with SE^ 1 ' = in the space of atomic 
orbitals. Then observe that, in view of the second relation in (|63|), this (A^— A^-dimensional 
space with SE^ = can be identified with the kernel of the matrix T^. One now finds from 
(|60|) that the (yet unknown) coefficients Uj^ should be associated with the solutions of 

5EfUj, =^V^ + E I V>S + (Wl<^j> (1) , (69) 

V V 

for \5<f)j)^ expressed in terms of Uj^ through \5<pj)^ = Y^^Uj^G^V^k^) (see fl67|)). Actu- 
ally, in (|69"D , it can be shown (using the property J2 V T^ v Uj v = 0) that the term r ^pOjl is 
irrelevant to this order, and therefore fl69|) is really an eigenvector equation for the vectors 
U Jfl : 

SEfUjp = ((K\V 2 \K) + (k,\V 2 G' 2 V 2 \k u )) U Ju . (70) 
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This equation may be used to determine the coefficients Uj^ and the energy shifts 5Ey>. As 
one can see from this consideration, our perturbative formalism can deal with essentially all 
situations regarding the bound-state problem with well-separated potentials. 

IV. CONCLUDING REMARKS 

In this paper we have presented a systematic perturbation theory for energy eigenstates 
when the potential of the system consists of two spatially well-separated pieces, under the 
assumption that complete energy eigenstates of the two local Hamiltonians are available 
for our use. Our perturbative development, an expansion in the number of wave-function 
stretching factors, is reminiscent of the multiple scattering series. Depending on whether 
the local Hamiltonians have (almost-) degenerate energy levels or not, different perturbation 
theories must be used. Especially, when the local Hamiltonians have degenerate energy 
levels, one obtains from our theory systematic higher-order correction terms beyond the 
predictions of the molecular orbital theory. The reasonably simple formulas we found for 
the leading correction terms, that is, fl2"0|) and (fn|) in the nondegenerate case and (|5l|)-(|5T|) 
(or flo7|)-(|T0|)) in the degenerate case, may have some immediate practical applications. 

Extension to the case with more than two spatially localized potentials (in fact even 
to the case of a lattice of potentials) should be straightforward. Also, if the degenerate 
atomic orbitals are present in association with certain symmetry in the system, one may 
utilize so-called symmetry-adapted linear combinations of atomic orbitals[l,2] to simplify the 
perturbation theory. But we have not made any systematic attempt in this direction. We also 
remark that if the local potentials happen to be not sufficiently well- localized (i.e., individual 
potentials have some long-range tails), certain rearrangements may become necessary with 
our perturbation series. This case deserves further study. One can also contemplate on a 
simple field-theoretic application: perturbation theory similar to the one given in this paper 
may be used to study the fermionic bound states associated with a soliton-antisoliton pair[0. 
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APPENDIX 



In this appendix we will first present the argument that shows why the terms in the right 
hand sides of ([18]) and fl38|) are expected to be of higher order than those terms in the left 
hand sides of the respective equations. We will then make estimates, by general argument 
and by considering explicitly the case of 5-function potentials, on how small the suppressed 
continuum contributions might be. 

Let us start with our equation (fL8D , used for the leading order approximation, and fl25|) for 
comparison's sake. The two candidates one wishes to identify as the leading-order expression 
of \5<j)k), i.e., G2V2IA;) according to flT8"D and G^V^/c) from Q25|), contain one wave-function 
stretching factor coming from the overlap of V2 and \k). Since 5Ek carries at least two wave- 
function stretching factors, the last terms in the right hand sides of flTBp and (^H) can safely 
be ignored in lowest order consideration. Here, the dangerous terms are the ones without 
SEk, i.e., 1^) in (JH) and i^G^V^lk) in (|25|) . Both have clearly one wave-function 

stretching factor from and one might expect that an additional suppression might 

result from the combination V1G2V2 or V2G1V2. The Green's operators G[ and G2 come 
with the sum over appropriate energy eigenstates (of the local Hamiltonians Hi, H 2 ), which 
include the continuum. As for the bound state contributions of one local Hamiltonian to the 
Green's operator, there should be such additional suppression (due to small overlap) if they 
get combined with the potential of the other local Hamiltonian. But, for the contribution 
to the Green's operator from continuum states which are not localized at all, one might 
not expect such suppression factor to show up, for these continuum states would apparently 
have more or less equal overlap regardless of the 'location' of the other potential. But this 
ignores the fact that one should really consider the net effect of entire continuum states. As 
will be discussed below, we get a very different picture after integrating over the continuum. 

Representing the continuous eigenstates of H2 by \q), we may express the continuum 
contribution of V1G2V2 as 

[^Mm . (71) 

Here, for well-localized potentials V\ and V2, the vectors Vi\q) = J dx\x)Vi(x)(x\q) and (q\V2 — 
J dy{q\y)V2(y)(y \ will receive nonnegligible contributions mainly from the regions around the 
respective potential centers, i.e., x = and y = L. Furthermore, outside the range of the 
potential V2, the function (x\q) may well be approximated by a plane wave. This implies 
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that Vi|g)(g|V2 comes with a phase factor 0)|Vi|g)(g|V2|j/(~ L)) ~ e l ^ L , which causes 
a destructive interference if the separation L is sufficiently large. Hence the term in (|7TD 
comes with desired additional suppression. On the other hand, an analogous consideration 
with the continuum contribution of V^G'^Vi does not lead to such a fast oscillating factor 
and so no suppression after summing over all corresponding continuum states. This explains 
why, for our leading order analysis, we can utilize flT8|), but not fl25|). 

By same reasoning as above, we expect that the continuum contributions for, say, (&IG2V2 
or Vi{G' x ) n V\ (with n > 2) be also suppressed. Suppression in the former case follows since 
Gi appears between the bound state \k), which is localized around the center of Vi, and 
the potential V2 (localized around x = L). As for the latter, the continuum contributions 
from {G'i) n give rise to a rapidly oscillating phase if sandwiched between V\ and V2. This 
suggests also a useful guideline in our consideration of higher order perturbation terms: to 
have the continuum contributions from G[ and G2 suppressed as much as possible, we had 
better iterate the relevant equation so that such Green's operator may take its place between 
V\ (or \k)) and V-j. In fact, we followed this guideline to obtain the expression for \5<f))( 2 \ 



i.e., when we proceeded from ([Tj|) to (]38f) . Of course, to confirm that (]38|) leads to the 
identification fl40|) for 150)^, we need to pay more careful attention to the order of various 
terms appearing in fl38|), and especially demonstrate the relative higher-order nature for the 
expression on its right hand side. For this, see below. 

First, with (|38|), look at the terms appearing inside the curly brackets on its left hand side. 



Based on (^), we know that \5<p k ) {1) is of order e^V 2 ™^! L ; the suppression factor here 
originates from the exponential tail of the bound state \k). With the next term G' 1 ViG2V2\k), 
we note that (aside from V2I&), itself of order e - RV 2m ' £fe L ) there is another small factor, 
say k,2, coming from V1G2V2. That is, this term is of order K 2 e _ «V 2m ' £fc L . The last term, 
5E k {G 2 ) 2 V2\k), is of order { e -b/ 2m \ e *\ L } 3 . Between G' 1 V 1 G 2 V2\k) and 5£ fe (G 2 ) 2 t> 2 |A;) we can 
not say generally which one is larger, because the magnitude of k 2 depends on the specific 
problem under study. Therefore, it is appropriate to include both terms in our second order 
approximation, under the understanding that only one term may well be dominant over the 
other in a given specific problem. Similar analysis can also be made for various terms in the 
right hand side of (|38|). If we denote the small factors emerging from V^G^Vi and Vi(G2) 2 V2 
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by K\ and k 2 , respectively, we here find 



V 2 G' 1 V l G 2 V 2 \k) ~ KiK 2 e~hPM^\ L , 5E k G' 1 V 1 G 2 V 2 \k) ~ /^(e^V 2 ^ L)3 



^t4(G 2 ) 2 t> 2 |A;) ~ 4(e"*V^i £)3 j (^) 2 (G 2 ) 2 F 2 |&) ~ (e^V 2 ^! L f. (72) 

This shows that the terms in the right hand side of (|38| ) are relatively of higher order, as 
compared with those terms appearing inside the curly brackets on its left hand side. [Note 
that this is true without any extra assumption about the relative ratios between K\, k 2 and 
k' 2 }. Hence our formula for the second order correction \5(pk)^ in (0) follows. 

In the above discussion, various small factors like Ki, k 2 have been introduced. How small 
are the continuum contributions associated with these factors? Consider V 2 G[Vi, with the 
related continuum contribution given by 

dxdy I dqV 2 {x) — ^—V^y), (73) 



2m 



where ip q {x) = (x\q), a continuous eigenstate of Hi. Here, for well-localized local potentials, 
the integral will get most of its contribution from the neighborhood of x, y ~ 0. For fixed x 
in the neighborhood, we may then extract the leading L-dependence from ip q (x + L) as 

ijj q {x) w A{q)e iA(k) e iq{L+x \ (74) 

where A(q), A(k) represent the amplitude and phase shift, respectively. With the form ( |T4"D 
used in (j73|), and after some careful study of the g-dependence including that from ip*(y), 
we notice that the given amplitude takes the form 

AkL 



dqf(q,x,yy kL , (75) 

where / corresponds to some regular function in q. For L very large, the order of magnitude 
for this integral can be deduced with the help of the Riemann-Lesbegue lemma[6]: if the 
n-th derivative of f(q) satisfies the so-called Dirichlet condition, then 

/_>)e-~0( I l T ). (76) 

Thus, for regular /, the integral should be smaller than any power of \. This strongly 
suggests that the typical large-L behavior of the integral (173) is that of an exponential 
suppression (i.e., vanishes like e~ aL , a being some positive constant). Remaining integrations 
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with respect to the variables x and y will not change this order estimate in any significant 
way, and so K\ is exponentially suppressed for large L. By analogous arguments one may 
demonstrate that k 2 and k' 2 , for large L, are also exponentially small. 

It is possible to give more precise large- L dependences for the factors Ki, K 2 and n 2 if 
a concrete problem is considered. For instance, we can compute these factors explicitlly 
when the problem is that of a pair of ^-function potentials as given in (^). Then, using the 
corresponding continuum wave-functions (see fl3lD ), we obtain 



2m 7 i7 2 r°° (x = L\q^)(q^\y = 0) /m 7l L 1\ _^nL f . 

Kl = 1 d « g2 + (z^)2 = ^ {-&- ~ 2 ) 6 ^' (77) 

K2 = -J^-j-^e-^-e-^l (78) 



7i-72 1 7i + 72 

,2. 



rnjl 



L 



, ^lil2 J 72 _vi22k 1 /WyiL 7i . ^_ . , 

™ I (7i -7 2 ) 2 7i(7i-72) V h 7i-72/ J 

As anticipated, we see the exponential dependences on L for these factors. We also observe 

from the results (|77|)-([79|) that, depending on the relative magnitudes of 71 and 72, it may 

be just one term that dominates the respective expression. Furthermore, as the result (|78|) 

for k 2 is used in our formula (fE]), we notice that the second term (proportional to SE^) is 

negligible compared to the first. Hence, for this example, we are allowed to write 



|c^ fc >( 2 ) = G'&G^k) = Kj^e-^G'^x = 0). (80) 
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